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Abstract 

The contribution of mesonic exchange currents to nuclear Compton 
scattering is investigated within the framework of a Fermi gas model of 
nuclear matter in the non-relativistic limit. The additional interaction 
between the nucleons is accounted for by including two- and three- 
body diagrams. As a test of this model, the enhancement constant k 
is calculated. The full correlators for the central and tensor part of the 
nucleon-nucleon interaction due to pion exchange are obtained and the 
energy dependence of the amplitude is investigated. The contribution 
of the A-excitation to the mesonic part of the Compton amplitude is 
calculated explicitely using an effective Hamiltonian in the static limit. 
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1 Introduction 



At present nuclear Compton scattering at intermediate energies is a highly 
advanced method to study nuclear and subnuclear degrees of freedom. Im- 
portant information can be obtained on the structure of the free nucleon and 
possible modifications of its properties due to binding. Compton scattering 
has been the subject of several recent theoretical and experimental articles 
|, |, |, |, I . One of the most promising subjects, which can be investigated 
in this way, is the contribution of mesonic exchange currents to the Comp- 
ton amplitude. This contribution is most important in the region between 
the energies of collective nuclear states and the A-excitation of the nucleon, 
overlapping with both at low and high energies, respectively. Experimen- 
tally it is a tedious task to separate these contributions, which entails that 
already for the analysis of the data important theoretical techniques are ap- 
plied, such as sum rules and dispersion relations (see e.g. |§, ^ 0)- Oii the 
other hand, the theoretical description involves numerous different meth- 
ods, depending on the energy region and the choice of nucleus (see Q and 
references therein). 

The experimental situation has improved significantly in the last years. The 
accuracy of Compton scattering experiments is now good enough to reliably 
extract some contributions to the scattering amplitude [^, |l^, |l^ . 

As recently discussed in Q the constant that reflects the strength of mesonic 
effects, the so-called "enhancement constant" k, enters into the description 
of nuclear Compton scattering at two stages. First, it describes the over- 
fulfillment of the dipole sum rule due to momentum-dependent potentials. 
Second, it is used to parameterize the additional term in the non-resonant 
amplitude due to Compton scattering by correlated proton-neutron pairs 
[p!3| . Thus, in addition to defining the low-energy behaviour of the Comp- 
ton amplitude, k is also present in the expression for the integrated pho- 
toabsorption cross section. The comparison of experimental data with the 
predictions on the grounds of dipole and quadrupole sum rules shows that 
K here contributes approximately 30 per cent |^] . 

In the present article we restrict ourselves to the consideration of heavy nu- 
clei. We assume that for large A'^ and Z the quasiclassical approximation 
can be applied, with A^ and Z being the numbers of neutrons and protons, 
respectively. We start from the Fermi gas model calculating the corrections 
to the corresponding nuclear wave function using perturbation theory with 
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respect to vr- and /9-meson exchange. These corrections are important, be- 
cause they determine an essential part of the nuclear correlation functions, 
giving a significant contribution to the mesonic part of the Compton am- 
plitude. In the course of this calculation it was realized that three-body 
corrections are as important as two-body corrections. 

Using the wave function obtained by this method, we constructed the 
mesonic contribution to the Compton amplitude and studied its dependence 
on the photon energy and the scattering angle. We also considered the terms 
arising from mesonic currents, when a A-isobar excitation is allowed as an 
intermediate state. 

The applicability of our results is restricted to photon energies less than 
approximately 150 MeV. This is due to our use of the meson-nucleon inter- 
action in the static limit (nonrelativistic approximation). 



2 Compton amplitude and Correlator 



In [ 14 1 it was suggested to write the total nuclear Compton amplitude Ta as 
a sum of three contributions provided by different physical mechanisms (see 
also [|l^] and |]l6|): the collective nuclear excitations (Giant Resonances), 
Tgr{u),6), the scattering by quasi-deuteron clusters, TQ£){io,6), and the nu- 
cleon Compton scattering, Tfyf{uj,9), where to is the photon energy and 9 
is the scattering angle. Here we will only deal with the first two of these. 
When the photon energy oj is small enough (less than pion mass), it is pos- 
sible to extract from these two contributions their resonance parts, which 
have non-zero imaginary parts at 9 = 0. Due to the optical theorem these 
imaginary parts are connected with the photoabsorption cross section. With 
the use of a dispersion relation one can obtain the real part to these reso- 
nance contributions in forward direction. Phenomenological arguments have 
been used to extend this resonance contributions to non-forward directions 
(see e.g. |l^]). Also a phenomenological description has been suggested 
for the non-resonant contributions (seagull parts) (see e.g. [|l4|). Unfortu- 
nately recent experimental results show discrepancies between this approach 
and the experimental data |[l|, ^. The strong restrictions on the resonance 
parts due to dispersion theory make them less model dependent than the 
non-resonant parts. Thus, in order to achieve an agreement with the exper- 
iments, we investigate the non-resonant parts with the use of microscopic 
considerations. 
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Due to the low-energy theorem, the Compton scattering amphtude at w = 
is equal to its Thomson limit 

TAiO,e) =-81-82^. (1) 

Here 8i and e*2 are the polarisation vectors of the incoming and outgoing 
photons, respectively, M is the nucleon mass and e is the proton charge, 
= 1/137. The low-energy limit of the resonance part of the Compton 
amplitude can be represented as 

R{0,e)= 81-82 — -{1 + ^), (2) 

where k is the so-called enhancement constant. This form (|2|) is related to 
the dipole sum rule for the photoabsorption cross section. 

The non-resonant part can be represented as a sum of two contributions: 
the kinetic seagull amplitude B{ijj,9) and the mesonic seagull amplitude 
Ai{uj,9). The kinetic seagull amplitude corresponds to the Thomson scat- 
tering by individual nucleons inside the nucleus: 

B{u;,e) = -8-,-8-2—Fp{A), (3) 

where 

F,{A) = j drp,{r)e^^'' (4) 

is the form factor belonging to the distribution Pp{r) of protons inside the 
nucleus and A is the momentum transfer, A = ^2 — ^i; with ki and k2 being 
the momenta of the incoming and outgoing photons, respectively. Therefore, 
in order to obtain eq. (||), it is necessary for the mesonic seagull amplitude 
to satisfy the relation 

ZN 

M{Q,Q) = -8i-82^e\. (5) 

Due to eq. ^ our calculation of k, as described in Section 3, can be compared 
with the value obtained from experimental photoabsorption data. This pro- 
vides a test of the low-energy properties of the amplitude in our approach. 

We consider mesonic contributions to Compton scattering in heavy nuclei 
and apply a Fermi gas model of nuclear matter for the nuclear wave func- 



tions. Such an approach was also used in Ref. pi, 1£] for the case of forward 
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direction. There it was shown that in such an approach all nuclear degrees 
of freedom can be taken into account by using a correlation function. Such 
a correlator consists of a central and a tensor part. In a pure Fermi model 
the tensor correlator is zero. However, the importance of tensor correlations 
was already pointed out in [^]. In its effect on the enhancement constant 
At, for example, the central correlator makes a comparatively small contribu- 
tion due to the strong compensation between the different diagrams. In [^] 
the results of numerical nuclear matter calculations for the correlators were 
used and compared to a similar calculation based on a variational principle 
pO[ . In order to obtain analytical results for the Compton amplitude and 
take into account not only the central part, but also the tensor part of the 
correlator, we calculated the corrections to the nuclear wave function, using 
the interaction between nucleons due to vr- and p-meson exchange as a per- 
turbation. It is clear that the accuracy of such an approach is limited. But 
we hope to reproduce qualitatively most of the features of the amplitude. 
Our analytical representation of the correlators allows us to study the de- 
pendence of the amplitude on different parameters (nuclear radius, N and 
Z), as well as its energy and angular distributions. 

We start our investigations from the static Hamiltonian for the pion-nucleon 
interaction (see e.g. pT|) 



H^NN = --(ff-V)(T-<j)), (6) 



m 



together with a minimal coupling to the photon field. In eq. (|^ m is the 
pion mass. For all coupling constants we use the same notation and values 



as given in [21|, in particular /^/(47r)=0.08. In eq. @ underlined symbols 
denote vectors in (cartesian) isospin space, while an arrow indicates a vector 
in coordinate space. 

The mesonic seagull amplitude A4 can be written in the following form: 



dQ 



The contribution T^*^^ of vr- meson exchange to Tj-,-, corresponding to the 
diagrams shown in Fig. (|^), is given by: 

('^) m? I Di Did2 D2d2 



4 



Q{Q\ ei • Qi £2 • Q2 _ Q2Q1 gi • g2 / j \ . 

where the following abbreviations have been used: 

2 ^ 
L'1,2 = (Q ± i?) + - Qi,2 = Q ± y , = g2 + 

with K = (ki + k2)/2. The correlator S*-' entering into eq. (|^ has the 
following general form: 

Here the summation with respect to a and h is performed over all nucleons, 
r^^^ = (n lb iT2)/2 are the isospin raising and lowering operators. The 
internal pion momenta were defined in such a way that photon and pion 
momenta add up to Q + A/2 in the left nucleon vertex and Q — A/2 in the 
right one. This gave us the possibility to write the amplitude Ai in the form 
(0), i.e. to extract the correlator (^). 

The correlator S"*-^ reveals some properties, which are interesting for investi- 
gation. Here we will study its dependence on the nuclear radius i?, on N and 
Z and on the momentum transfer A. We will do this for the simplest case, 
i.e within the framework of a pure Fermi gas model. Then the correlator 
only has a central part: S^^ = SqS^-^ , where 

Sc = -2j j dxidX2 e-*(-i+-=^)A/2 ^^{x.-S2)m-Q-p,)_ (10) 

The ranges of integration for the nucleon momenta pi and p2 are the proton 
and neutron Fermi spheres, with the radii 

#=2.27m (1)'^' , # = 2.27m (^y^' , A = N + Z. 

In eq. (jl^) the integrations with respect to xi are taken over a sphere with 
the radius R = 1.2 A^^^ fm. Straightforward integration gives 



"6^ 



00 



sc = — ^ E (2^ + m (x) 



1=0 
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dv v[u — 2) 



+2(2 -cf I duu'gf\u)g'^'{i.) 



fill 



with X = (Qi • Q2)/( Qi Qi), a = pf,b = p'p' and c = 2(a - 6)/(a + b). 
The functions (7j-'^ can be expressed via spherical Bessel functions ji, 



An) 



ji{vR)QiRji.i{QiR) - ji-i{vR)vRji{Q^R) 
{vRY - {Q:,RY 



(12) 



Let us investigate the dependence of the correlator on the nuclear radius R. 
For simphcity we consider the case A = and N=Z. Then we have 



Sc\ 



A=0. 
N=Z 



s. 



(/) 

c 



A {Ra) 
2" Svr 



l + -)x 



+1 

j dX F^{Ra^Jq^ + -2quX), 



(13) 



-1 



with q = Q/pF and pp = p^P = p^p^ = 1.8m. The function F is given by 



F{x) 



3 / sin X 

x'^ \ X 



cos X 



(14) 



The superscript (/) for the correlator denotes that this result corresponds 
to a finite nuclear radius. Note that in a Fermi gas Ra = {9ttA/8)^/^. If R 
tends to infinity, one obtains the asymptotic form of (p!^): 



s: 



(0) 

c 



A 
2" 



1 + 



e{2 



(15) 



which is a well-known result (see e.g. [22]) of the Fermi model. Here 0{x) 



is the step function. In Fig. (g), eq. (13) is compared with eq. (|l5|) for A=50 
and ^=200. One can see that even in the case of ^=50 the use of the 
asymptotic form for the correlator is justified, even though in the region of 
very small Q the curves differ significantly. 
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We now consider the correlator Sc from eq. (11), when A=0, R^oo , but 
Nj^Z. Then we have, for N > Z: 



a-Q) (16) 

e{a + b-Q)e{Q + a-b) |. 

Let us introduce the parameter x = {N — Z)/A. In Fig. (^) the quantity 
Sq^ is shown as a function of Q/pf for different values of x. One can see 

that even for large x ~ 0.5 the difference between and is essential 
only in the region of very small Q. In our further analytical calculations we 
will neglect this dependence and put N=Z=A/2. In addition to this we will 
consider the limit of large i?, but i?A ~ 1. Then it is easy to show that the 
correlator Sc can be represented in the following form: 

Sc = SP FiRA), (17) 

where F and S^^ were defined in eqs. (0) and ([T5[), respectively. For the 
corrections to the correlator, which will be investigated in the next section, 
the dependence on the momentum transfer A is also determined by the form 
factor F{RA). 



~(0) (g + by 

c 2 ' " 



1 

+- 



1 - 



(a3 + 63) 
a + b) 



Q2_3(^_fc)2 - 

Q{a + b) 



3 Corrections to the Correlator 

As was mentioned in Section 2, in order to obtain the correct behavior of 
the Compton amplitude, it is necessary to take into account not only the 
central part of the correlator, but also its tensor part. Therefore, in this 
section we calculate the correction to the wave function. This can be done 
with the use of usual perturbation theory up to second order. In this case we 
have to consider two-body and three-body effects with respect to correlated 
nucleons. Again all dependence of the correlator on the scattering angle is 
given by the form factor F(i?A), as pointed out in the last section. So, let 
us consider eq. ^ in the case of A=0, which without change of notation we 
will represent in the following form: 

S^i = Sc S'^ + St f^, (18) 
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where 



3Q^ 



(19) 



The diagrams, which correspond to the two-body part of the correlator 
correction, are shown in Fig. (Q). Exphcit evaluation of these diagrams leads 
to the following correction: 



q}3 



4:Mf^ fdpidp2 j iQ'Q^ ^ 2rW^-r^6'^ 



(27r)6 1q2+^2 



+ 



n{pi)n{p2) 1 - n{p2 + Q) 1 - n{p 



-Q) 



where 



Q 



PI-P2-Q 



(20) 



and V is the nuclear volume. The function n{p) = 9{pf — \p\) is the occupa- 
tion number in the Fermi gas model. Due to the symmetry of the integrand 
the product of four occupation numbers can be omitted and one can replace 
[1 - n{p2 + Q)][l - n{pi - Q)] by [1 - 2n{pi - Q)]. Eq. (|0|) contributes to 
both, the central and tensor part of the correlator. We can separate these 
two contributions with the help of the identity 

and the fact that due to the (spherically symmetric) integrations with re- 
spect to the nucleon momenta pi and p2 the expression (2rV-' — can 
be replaced by 

- f2<5^i . 

Analytical integration of eq. (|2^) with respect to pi and p2 leads to the ex- 
plicit formulae presented in the appendix. It is important to realize that, in 
order to fulfill the requirement of gauge invariance, one has to evaluate also 
the set of three-body diagrams displayed in Fig. (P). As a result all contri- 
butions of the second order in the meson-nucleon coupling are taken into 
account. In accordance with the Feynman rules one obtains the following 
general expression for the three-body contribution to the correlator: 




•^(2) 



l6Mf 



V 



hdp2 ) 2Q'Q^ 



(27r)6 [ Q2 + ^2 



+ 



2t^t-^ — j^'^^^j 
r ^ -|- 
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Q-f ' 

Here the product of four occupation numbers has also been omitted, as in 



eq. (2C). The result of the explicit analytical integration for eq. (|2l| ) can be 



found in the appendix. As before, the correction (21) contributes to the 
central part of the correlator, as well as to its tensor part. Thus, at this 
stage of our investigation the correlation function S^^ has the following form: 

S^j = 3^6'^ + S'^^^+Siiy (22) 

Note that central and tensor parts of the corrections can be extracted by 
using the definitions ( [l^ and (|l9|). In Fig. (^) the full central part of the 
correlator as a function of Q/pf is shown, together with the zeroth order 
approximation S^^ and the one, where only two-body correlations have been 
taken into account. It can be seen that almost over the whole range of Q the 
asymptotic form of the central correlator is modified, although the overall 
shape remains the same. 

Let us now pass on to the discussion of the tensor correlation and the differ- 
ent contributions to it. Fig. (0) shows the tensor correlator with and without 
three-body effects. Here a strong modification due to three-body correla- 
tions occurs at comparatively low Q, where we found a strong damping that 
is in fact essential for obtaining reasonable physical results. A good test of 
the model developed so far is the calculation of the enhancement constant 
K. With the use of eq. (^), together with eqs. (|^) and (^), we obtain 



oo 

16 M 



K = — — 




Sc-2s^St 4s2(5c + 25r) 



{s^ + iy 3(s2 + l) 



s'^ds, (23) 



where the functions Sc and St depend on Q = ms and contain the two- and 
three-body corrections, which are given in the appendix. Our calculation 
is restricted by the use of the nonrelativistic approximation. We, therefore, 
cannot apply our results at large Q/pf- One can see from Fig. (0) that the 
contribution of the high-Q region is not negligible. It is known that taking 
into account only the one-pion exchange leads to a wrong behavior in the 
tensor part of the nucleon-nucleon potential at small distances (i.e. at large 
Q) (see e.g. [^]). Thus, it is necessary to include /o- meson exchange in the 
calculation of the tensor part of the correlator. This can be done ||2^ by 
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replacing St by: 



St = St 



1 



(24) 



g2 + mj/pl, 

where nip is the /o-meson mass and the following relation for the pNN- 
coupling constant has been used: 



fp 



rrin 



-2 



Taking into account the p-meson improves essentially the behaviour of our 
correlator at Q ~ p_p and higher (cf. Fig. (^)). Unfortunately the uncertainty 
connected with the choice of the upper limit in the integral representation 
for K, eq. (p3|), still exists. The integrand for the s-integration in (23) is 
shown in Fig. (|9|). We estimate the accuracy of the result obtained for k 
to be of the order of 20 per cent. Introducing a cut-off at Q/pf = 2.7, 
we obtain the numerical values for the different contributions to k. Using 
correlator gives ko=0.2, in agreement with [^]. With the full 
central part of the correlator, one finds kc=0.15. The contribution of the 
tensor part is kt=0.86. So, the total value for the enhancement constant is 
K = Kc + kt = 1.01. 

Mind that all results for k, obtained in this section have been calculated for 
infinite nuclear matter with N=Z. In Fig. (|3|) the influence of the difference 
between and Z has been shown for the lowest-order form of the correlator. 
It is also interesting to study the enhancement constant k, which does not 
depend on the momentum Q, as a function of N and Z. We will do that 

for the simplest case of only the S'^'* -contribution to k (i.e. kq). Let us 

:{0) 



put in eq. (P3) St = 0, Sc = Sq (from eq. (|16D ) and multiply the rhs of 
( p3| ) by A/[2N) in accordance with the definition (^). After straightforward 
integration one obtains 



AM fV 



-ab +- (a^ + 771^ + b'^] In 



\ m J 



+ 1 



(25) 



where a = 
is shown. 



Pp'\ as before. In Fig. (|To|) the {N, Z)-dependence of kq 
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4 Energy and Angular Dependence 



In the previous section we used the full central and tensor parts of the nuclear 
correlator to calculate the enhancement constant k. Let us pass now to the 
consideration of the energy dependence of the amplitude (0) . We start from 
the low-energy limit of the amplitude and take into account terms of the 
order of uP' Ivr? . Making the expansion of T'^i^^y eq. (^), with respect to K 

and A up to second order, contracting with b^^ and t*-' and, finally, averaging 
with respect to the angles of Q, we obtain within our accuracy 



M = ^F{R/\) 



where 



4M 



a 



a uj 



£i • £2 H n ei X fei 



£2 X h 



(26) 



^^^^ ' s^ds 



-2St 



{Sc + 2St) 



lOs^ + 39s^ + 60^2 + 15 
15(1 + 



3s^ + 5 
5(l + s2)3 



and 



P 



4M/2 



s ds 



{Sc + 2St)^47T^^ + 2St ^ 



3(1 + ^2) 



(1 + S2)3 



(27) 



(28) 



The constants a and /? can be interpreted []18|, |14[ as the modification due to 
medium effects of the electromagnetic polarizabilities for the free nucleon, 
measured in units of (^=6"^ / {4Mm'^)=7 .5 x 10~^fm^. Including this coeffi- 
cient we obtain the following numerical results: a=(^a=— 3.23 x 10~^fm^ and 
/3=C/3=1.73xl0-4fm3. The use of only S^^U'^ leads to ao=-4.5x IQ-^fm^ 
and ^0=1-95 X 10~^fm^. The sum a-|-/3=— 1.5x 10~^fm^, which is obtained 
from the consideration of forward direction scattering, gives a correction of 
the order of 10 per cent to the corresponding value for the free nucleon. 
The polarizabilities contribute to the backward scattering amplitude in the 
combination /3— q=5x 10~^fm^. Here the relative correction is essentially 
larger than in the case of forward direction. 



Note that the expression in brackets in eq. ( pq ) does not contain an electric 
quadrupole contribution. Some E2 strength occurs due to the expansion of 
the form factor F{RA). However, this contribution explicitely cancels with 
its resonance counterpart, in accordance with the low-energy theorem. 
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It is interesting to compare the exact energy dependence of our amplitude 
with its low-energy form (^). In Fig. ([Tl|) this dependence is shown for 
the cases of forward and backward scattering. One can see that up to an 
energy of 100 MeV all energy dependence is reproduced by eq. (^). We 
can, therefore, apply the low-energy form of the amplitude at all scattering 
angles. 



5 Virtual A-excitation 



Let us consider now the contribution of a virtual A-isobar excitation to the 
mesonic part of the Compton amplitude. General properties of such an effect 
were already discussed in |24]. Here we have the possibility to consider this 
contribution in more detail with the use of the correlation function obtained 
in Section 3. In the static limit the corresponding Hamiltonians, which 
determine the interaction, are of the form (see e.g. pT|) 



H, 



and 



■yNA 



m 



V X A] T. 



m — 



(29) 



(30) 



with the hermitian conjugate to be added in both cases. Here S and T are the 
l/2-to-3/2 transition operators in spin space and isospin space, respectively. 
Evaluating the diagrams given in Fig. ( [T^ ) on this basis, one obtains the 
following expression for the isobar contribution to the tensor T^^ entering 
into eq. (^: 



7] 



di 



£*! X {k2 X £2) - £2 X (^1 X £1) 



(31) 



with 



ha 



2{Qa ■ £a)Qa 



da 



la = {Q + K) X {ka X Ea) 



and the coefficient 



SQe'^fAf-yNAf 



In (|TD the same abbreviations have been used as in eq. (g) and, in addition, 
0, = Ma — M is the mass difference between the A-isobar and the nucleon. 
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Also, the symmetry with respect to the substitution Q —Q, which is due 
to the integration in eq. has been used to bring T^'^-^ into the form 

(|3l|). As before, we consider the low-energy limit of the isobar contribution 



to the amplitude, which again can be represented in the form (2(:). In this 
procedure we find aA=0 and 

. 16M/aA7Va/ f 2 , '^s^St + {s^ - 3)Sc 

= 27.-n J ' WTTy • ^''^ 

Inserting the numerical values for the coupling constants /a =2/ and 
/7ArA=0.35, we obtain /3a=0.09 and, correspondingly, for the contribution 
of the A-isobar to the nucleon magnetic polarizability in the usual units 
^^=0.7xl0-'^fm3. As before, the energy dependence of the amplitude Ma 
discussed in this section is well described by its low-energy limit. At energies 
up to the pion mass Ma is less than 0.1 in units of Ae'^ /{AM). 



6 Conclusion 



In this paper we have tried to understand qualitatively the influence of 
different correlations inside heavy nuclei on the mesonic contribution to 
Compton scattering. Our model approach gives us the possibility to perform 
a large part of the calculation analytically and, as a result, to understand in 
more detail the interplay between the different physical contributions as well 
as the influence of various parameters. We started from a simple Fermi gas 
model and then investigated the correlations, which appear in such a gas due 
to vr- and p-meson exchange between the nucleons. In the calculation of these 
corrections, three-body diagrams have been found to give a contribution of 
the same order of magnitude as two-body diagrams. For example, strong 
compensations between two- and three-body contributions to the correlator 
occur in the region of relatively small pion momenta. 

We obtained that the tensor part of the correlation function has as important 
an effect on the Compton amplitude as its central part. In particular it 
increases essentially the enhancement constant k. Such nuclear correlations 
have also been found to be significant in the modification of the nucleon 
electromagnetic polarizabilities via mesonic currents. In addition we found 
that it is not unnatural to introduce, in the analysis of experimental data, 
a rather high modification of (3— a. 
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The contribution of the virtual A-isobar to the magnetic polarizabihty has 
been found to be less important than that of the purely pionic amplitude. 
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Appendix: Analytical expressions for the central 
and tensor parts of the correlator 

Here we present the explicit expressions for the correlator corrections, which 
have been used to obtain the figures shown in Section 3. The two-body cor- 
rection Slis , eq. (|20|), leads to the following additional terms in the correlator: 




-C [Gieiq-2)+G29{2-q)] 



(33) 



for the central part and 




for the tensor part of the correlator. In (33) and ( |3^ we used the following 



abbreviations: 



The functions Gi and G2 are of the form 



1 = h [n\2 — q 

60 120 15q 12 24 480 ' ^ 
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+ 



and 



Go 



12 



29 
240' 

2 



240 



ln(g) + 



1 

15^ 



12 



24 480 



ln(g + 2) (35) 



15g 6 



320 



+ 



1 

15g 



ln(2) + 

The integrand w{q,x) is given by 



1_ g2 

16 96 

1 1 

TSg 16 ~ 



1280 



96 



1280 



ln(2 - q) 

ln{2 + q). (36) 



w{q,x) 



1 

32^ 



1 



2x9 



x(x^ + 1) 



2'^ 



2 -1)2 In 



x + l 



X 



1 



(37) 



The external factor is C = {3M ppp)/ {n^m'^). 

It is clear that the main difference between the calculation of two- and 
three-body correlations occurs in the structure of the product of occupation 
numbers n{p) (cf. eqs. (p0|) and (^)) It is convenient to perform some 
of the integrations in ( ^l]) by representing the integrand as the one for the 
two-body contributions plus some additional terms. This procedure leads 
to the following expression for the correction 5^2) due to three-body graphs: 



5g) = 2 5[l) + 55^^ 



2CA q^ (G2 - G 



(1) 



f + 72 



+ 6'^) e{2-q), (38) 



with 




^-1 i6-+f^) 



In 



1 + 



(g-Hx)2 -K72 

(g — x)2 -|- 72 



72 {6'^ + 1'^) 
6q (x2 -|- 7^ — g'^ 



+ 



2g2 



72 -I- q'^x'^ 
Aqx 



In 



(g -f x)2 72 

(g — x)2 -|- 72 



(39) 



15 



Here 7 = m/pp. The coefficients of J*-' and t*-' in eq. ( p8[ ) are the functions 
Sq^ and S*-^^^ used in Section 3. 
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Figure Captions 



Figure 1 Typical diagrams contributing to T^^y The wavy lines denote 

photons and dashed lines denote pions. The amputation indicates that T^^ 
contains only the nucleon vertices, but not its wave functions. 

( f) 

Figure 2 Comparison of the exact correlator S*^ with its asymptotic form 

Sq'^ obtained for i? — > oo (curve a)). In S(j^ A is set to 50 (curve b)) and 
200 (curve c)). 

Figure 3 Dependence of S^^ on Q/pf for different x = [N — Z) / A in units 
of NZ/A. The values used for x are a) x =0.5, b) x =0.25, c) x =0.1 and 
d) X =0. 

Figure 4 Set of diagrams, from which the two-body correction to the cor- 
relator is extracted. Here the same symbolic abbreviation is used as in Fig. 
The nucleon spin projections are denoted by Aj. 

Figure 5 Three-body diagrams yielding an additional correction to the cor- 
relator, which is of the same order in f /m as the two-body diagrams shown 
in Fig. The notations are the same as in Fig. ^) and Fig. 

Figure 6 Central part of the correlator a) with three-body corrections, b) 
without three-body corrections and c) Fermi correlator without any correc- 
tions. 

Figure 7 Tensor part of the correlator a) with three-body corrections and 
b) without three-body corrections. 

Figure 8 Comparison of the tensor part of the correlator a) with and b) 
without the p-meson contribution. 
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Figure 9 Integrand for the calculation of k shown as a function of Q/m: 
a) the full integrand, b) contribution of the central part and c) contribution 
of the tensor part of the correlator. 



Figure 10 Dependence of kq on N and Z. 



Figure 11 Energy dependence of the amplitude Mn- The exact dependence 
(full line) is compared with the low- energy form (dashed line) for scattering 
a) in forward direction and b) in backward direction. A^^r is shown in units 
of-Ae^/{4M). 

Figure 12 Set of diagrams, which arises from allowing a A-isobar excita- 
tion as an intermediate state in the diagrams displayed in Fig. 
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